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Abstract We study quasi-geostrophic turbulence and plasma drift turbulence within 
the Charney-Hasegawa-Mima (CHM) model. We focus, theoretically and using nu- 
merical simulations, on conservation of zonostmphy and on its role in the formation 
of the zonal jets. The zonostrophy invariant was first predicted in (TKH in two spe- 
cial cases - large-scale turbulence and anisotropic turbulence. Papers fll, '21 also 
predicted that the three invariants, energy, enstrophy and zonostrophy, will cascade 
anisotropically into non-intersecting sectors in the A:-space, so that the energy cas- 
cade is "pushed" into the large-scale zonal scales. In the present paper, we consider 
the scales much less than the Rossby deformation radius and generalise the Fj0rtoft 
argument of O] |2| to find the directions of the three cascades in this case. For the 
first time, we demonstrate numerically that zonostrophy is well conserved by the 
CHM model, and that the energy, enstrophy and zonostrophy cascade as prescribed 
by the Fj0rtoft argument if the nonlinearity is sufficiently weak. Moreover, numer- 
ically we observe that zonostrophy is conserved surprisingly well at late times and 
the triple-cascade picture is rather accurate even if the initial nonlinearity is strong. 



1 Introduction and the model 

Zonal jets are prominent features in geophysical fluids, e.g. atmospheres of Jupiter 
and Saturn |l4l|5]|3 and the earths' atmosphere fl0l [3l and oceans |I3]|6]|3. Zonal 
jets have also been observed in fusion plasmas fTSl. Zonal jets are important be- 
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cause they can suppress the small-scale turbulence and block the transport in both 
geophysical settings ifTTI and in plasmas lfT2l [l4l [TSl . 

Two main zonal jet generation mechanisms considered in the literature are the 
modulational instability |16,8, 17,9 18, 19l|20l|2T| and the anisotropic inverse cas- 
cade ll22ll26ll23lfT2l . The inverse cascade mechanism brings the energy from initial 
small-scale turbulence to the large-scale zonal flows in a step-by-step (local in the 
scale space) transfer mechanism similar to the inverse cascade in 2D Navier-Stokes 
turbulence ||241 1251 . On the other hand, the modulational instability brings energy 
to the large zonal scales directly, skipping the intermediate scales. It appears to be 
more relevant if the small-scale turbulence has a narrow initial spectrum, whereas 
the cascade picture is more accurate when the initial spectrum is broad. We will 
report on our recent analytical and numerical studies of the modulational instability 
elsewhere 121], whereas the present paper will focus on the turbulent cascades. 

The mechanism for an inverse cascade in geophysical quasi-geostrophic (QG) 
turbulence and in plasma drift turbulence is quite similar to the one of the 2D Navier- 
Stokes turbulence l f24l l25l . but the presence of the beta-effect (non-uniform rota- 
tion of geophysical fluids or plasma inhomogeneity) makes this cascade anisotropic 
leading to condensation into large-scale zonal flows rather than round vortices 
|22, 26, 23 1 . 

In this paper we will follow the approach put forward in papers 1 1 2 1 which is 
most relevant (and asymptotically rigorous) when the QG/drift turbulence is weak. 
In this case the turbulence is dominated by waves which are involved in triad inter- 
actions. These three-wave interactions are shown to conserve an additional positive 
quadratic invariant. This invariant and the other two quadratic invariants, the en- 
ergy and the potential enstrophy, are involved in a triple cascade process, which can 
be described via an argument similar to the standard Fj0rtoft argument originally 
developed for the 2D Navier-Stokes turbulence Ii24il . It was found that each of the 
invariants is forced by the other two to cascade into its own anisotropic sector of 
scales and, in particular, the energy is forced to cascade to long zonal scales. Con- 
sidering its important role in the zonation process, hereafter we will label the extra 
invariant as zonostrophy. 

On the other hand, work ||T1[2 was limited to considering either very large scales 
(longer than the Rossby deformation radius or the Larmor radius) or to the scales 
which are already anisotropic and are close to zonal. Besides, the conservation of 
the extra invariant is based on the weakness of nonlinearity and on the randomness 
of phases (conditions of validity of the wave kinetic equation), which, even when 
present initially, can break down later on during the zonation process. Thus, numer- 
ical checks of robustness of the zonostrophy conservation were needed, and they 
have not been done until the work reported in the present paper 

Soon after papers the zonostrophy invariant was generalized to the whole 

of the A:-space in 1271 . This was a significant achievement because the extra in- 
variant of such a kind appeared to be unique for Rossby/drift systems and is not 
observed in any other known nonlinear wave model |28 1. Besides, its conservation 
has revealed interesting geometrical properties of the wave dispersion relation l(28l . 
Unfortunately, the general expression for zonostrophy appeared to have a form for 
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which the Fj0rtoft argument cannot be used (not scale invariant, not sign-definite). 
However, an alternative zonation argument was put forward in |29|. 

In the present paper we will focus on the special case when the scales are much 
smaller than the deformation (Larmor) radius, which is the most important and fre- 
quently considered limit (at least in the GFD context). Taking the respective limit 
in the general zonostrophy expression obtained in (TT\, we obtain the zonostrophy 
expression for such small-scale turbulence and show that it is positive and scale 
invariant. The latter observation is extremely important because it means that we 
can once again apply the generalized Fj0rtoft argument developed in |[Tl|2l, which 
is done in the present paper. Note that the Fj0rtoft argument of ||Tl|2l is somewhat 
preferential over the argument presented in Il29l because it predicts not only zona- 
tion but also the anisotropic A:-space flow paths of the three invariants during the 
zonation process. 

Having obtained these analytical predictions, we then proceed to direct numeri- 
cal simulations (DNS) of the QG/drift turbulence to (i) test the conservation of the 
zonostrophy for different levels of initial nonlinearity, and (ii) test predictions of 
the generalized Fj0rtoft argument by tracking in time the transfer path of the three 
invariants in the fc-space. As a result, we confirm conservation of zonostrophy and 
its important role in dkecting the energy to the zonal jet scales. 



2 Charney-Hasegawa-Mima model 

The reason why geophysical quasi-geostrophic (QG) flows and plasma drift tur- 
bulence are often mentioned together (in particular when discussing the zonal jet 
formation) is that some basic linear and nonlinear properties in these two systems 
can be described by the same PDF, the Charney-Hasegawa-Mima (CHM) equation 



where yf is the streamfunction, /3 is a constant proportional to the gradient of the 
horizontal rotation frequency or of the plasma density in the GFD and plasma con- 
texts respectively. In the GFD context, the jc-axis is in the west-east and the y-axis 
is along the south-north directions respectively. In plasmas, the y-axis is along the 
plasma density gradient and the x-axis, of course, is transverse to this direction. 
Here, we consider only small-scale turbulence with scales much smaller than the 
Rossby deformation radius in GFD and the ion Larmor radius in plasma contexts 
(this has already been taken into account in the CHM model, ([T]i). 

Let us put our system in a periodic square box with side length L and introduce 
the Fourier transform of the streamfunction. 



1301 EH: 



d,A V/ + jS^.i-i/ + {dx\l/)dyA Xj/ - {dy\l/)d_^A 1/ = 0, 



(1) 
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where k — {k^^ky) is a wavevector in the 2D plane. In ^-space the CHM equation 
([l]) becomes: 

^rV>k = ;(SkV'k+:^ ^ r(k,ki,k2)v/k, V^kj, (2) 

^ ki+k2=k 

where 

C0k = -§^ (3) 

is the dispersion relation for the frequency of linear waves (Rossby or drift waves in 
the GFD and plasma contexts respectively), k= |k|, and 

(ki Xk2),(/fc?-/t?) 

r(k, ki , kz) = - ^ (4) 

is a nonlinear interaction coefficient. 



3 Conservation of energy and enstrophy 

It is well known that the CHM equation, ([T]) conserves the energy and the enstrophy 
which in physical space are defined respectively as 

^=\\ (Vv/)' (5) 

and 

Q = yiAxi/fdx. (6) 

In Fourier space these conserved quantities can be expressed in terms of the wave 
action 

E = j \(Ok\ni,dk (7) 
^ = J kxtik dk. (8) 



We have 



and 



4 Conservation of zonostrophy. 



The energy and the enstrophy are exact invariants of the small-scale CHM model 
([T]). Zonostrophy is an exact invariant of the kinetic equation: 
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«k = J |ymr5(ki+k2-k)5(«(ki) + a)(k2)-a)(k))x 

[n(ki)n(k2) -2n(k)n(ki) sign-n(k)«(k2) sign (^;c^2a')] dkidk2, 

where 

Thus, it should be clear that the zonostrophy is only proven to be a conserved quan- 
tity under the same conditions for which the kinetic equation (|9]l is valid, namely 
weak nonlinearity and random phases. It is presently unknown if the zonostrophy 
conservation extends to a broader range of situations or not. Thus, the numerical 
tests of zonostrophy conservation are crucial, and this is one of the aims of the 
present work. 

In terms of the wave action spectrum, zonostrophy Z can be written as 



j Ck«k 



dK (9) 



where function i^t is the density of zonostrophy in the A;-space which satisfies the 
triad resonance condition 

C(k) = C(ki) + C(k2) (10) 

for all wavevectors k,ki and k2 which lie on the resonant surface given by the 
solutions of the wavevector and frequency resonance conditions, 

k = ki+k2 (11) 

and 

C0(k) = C0(ki) + C0(k2). (12) 

Expressions for i^k were first found in 1 1 , 2] in the special cases of large-scale tur- 
bulence {pk <C 1 where p is the Rossby deformation radius in GFD or ion Larmor 
radius in plasmas) and for anisotropic turbulence ( A:, <C ^y)^ ^fter which a general 
expression was found for all k's in 1271 . This expression is 

C(k) = arctan^ ^ arctan^ . (13) 

pk^ pk^ 

This expression can change sign and it is not a scale-invariant function of the 
wavevector components, which means that one cannot use the generalized Fj0rtoft 
argument of 1 1 , 2 1 to find the cascade directions in such a general case. 

To find zonostrophy in the special case of small-scale turbulence (which is con- 
sidered in the present paper) we have to take the limit pk-^ °o in the general expres- 
sion ( 13 1. It turns out that naively taking limit leads to an already known integral. 



the energy and one has to go to further orders in the Taylor expansion of (1 3 1 until 
we reach an expression which is independent of E and £2. Interestingly, such an 
independent invariant appears only in the fifth order, and the derivation details are 
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given in Appendix A. The result is 

C = - jim 0^ (C - 2V3CO/PP) = |J ie + 5k/). (14) 

The integral (j9]l with the density ( [T4| is an exact invariant of the kinetic equation (j9]l 
and thus it is an approximate invariant of the small-scale CHM equation ([T]l. Later 
on we will examine numerically the precision with which this invariant is conserved. 



Expression ( 14 1 for the new invariant in the small-scale limit, pk^ \, allows us 
to explicitly see that the invariant's density is strictly positive and scale-invariant, 
which is good news because once again one can use the generalized Fj0rtoft argu- 
ment of to find the cascade directions of the three invariants: the energy, the 
enstrophy and the zonostrophy. 



5 Triple cascade behaviour. 

5.1 Dual cascades in 2D Navier-Stokes turbulence. 

Before we present the generalized Fj0rtoft argument, we will remind the reader of 
the classical Fj0rtoft argument leading to the prediction of the dual cascade be- 
haviour in 2D Navier-Stokes turbulence ll24ll . This will be instructive because by 
now there exists several versions of such an argument and the one we use here is not 
necessarily the most famihar. 



Energy Enstrophy 




K- Kf K+ 

Fig. 1 Dual cascade behaviour in 2D Navier-Stokes turbulence. 



Let us consider 2D turbulence described by the Navier-Stokes equations excited 
at some forcing scale ^ and dissipated at very large (~ A:_ ^ fco) and at very 
small scales k+ ^ ko), see Fig.[T] The conservative ranges between the forcing 
scale and the dissipation scales, k_ ^k<^ko and ko <^ k <^ k^ sae called the inertial 
ranges. The two conserved quantities in this case, in the absence of forcing and dis- 
sipation, are the energy and the enstrophy which are given by the same expressions 
as for the small-scale CHM model, namely by (jSj) and (j6]) in the jc-space and by 
^ and (|8]l in the A:-space. In the presence of forcing and dissipation in steady-state 
turbulence, the rate of production of the energy and the enstrophy by forcing must 
be exactly the same as the dissipation rates, e and rj respectively. Our task now is to 
determined where exactly E and £2 are dissipated, at ~ k- or at ~ 
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First, we note that the A:-space densities of E and Q, differ by a factor of k^, and 
therefore the energy dissipation rate £ is related to the enstrophy dissipation rate 77 
as 7] ko^e. Now, let us suppose ad absurdum that at k+ the energy is dissipated 
at a rate comparable to the rate of production at the forcing scales, i.e. ^ e. But 
this would mean that the enstrophy would be dissipated at the rate ^ k^e which is 
impossible in steady state since this amount greatly exceeds the rate of the enstrophy 
production rj ko'E. Thus we conclude that most of the energy must be dissipated 
at the scales ~ i.e. that the energy cascade is inverse (with respect to its direction 
in 3D turbulence). 

Similarly, assuming ad absurdum that most of the enstrophy is dissipated at ^ A:_ 
would also lead to a conclusion that the amount of energy dissipated is much greater 
than the energy produced, which is impossible. Therefore, the enstrophy cascade is 
direct, i.e. it is dissipated at wavenumbers ~ k^ which are greater than the forcing 
scale ko. 



5.2 Triple cascades in CHM turbulence. 

From the Fj0rtoft argument presented above, the reader should notice that the quan- 
tities that determine the cascade directions are the A:-space densities of the invariants 
or, more precisely, the scaling of the ratios of these densities with k. It is also im- 
portant that these fc-space densities are positive (as is the case for E,£2 and Z in 
our situation, see (j7]i, ([8]l and ( 14 1) because otherwise a large positive and a large 



negative amount of the same invariant could be produced in different comers of the 
A:-space, with the net result zero and the Fj0rtoft argument would not work. Further, 
it is important that all of the invariants have the same scaling with respect to the 
turbulence intensity so that their ratios are functions of k and not of the turbulent 
intensity. This condition is satisfied, for both Navier-Stokes and in the CHM model 
since E,Q and Z are linear in ni^, see (|7]i, ([8]) and ( [T4| . 

The ratio of the A;-space densities of Q and E isP which allows us to conclude 
that Q must go to ^ ^ ko and E must cascade to k ^ ko. In other words, each of 
the invariants cascades to the scales where its density is dominant over the density 
of the other invariant. This argument is easily generalizable to the CHM model and 
the respective three invariants, E,Q and Z. Now we have three invariants, and the 
cascade picture would necessarily be anisotropic (it is impossible to divide the 2D 
k-space into three non-intersecting cascade regions in an isotropic way). 

Let us suppose that turbulence is produced near ko = {kox,koy) and it can be 
dissipated only in regions which are separated in scales from the forcing scale, i.e. 
either at k^ ko (short scales), or at k^ <C ko^ (nearly zonal scales), or at ky ^ koy 
(nearly meridional scales), see Fig. (j2]l. Then each of the invariants (e.g. E), must 
cascade to the scales where its density is dominant over the densities of the other 
two invariants (e.g. Q and Z). The boundaries between the cascading ranges lie on 
the curves in the A:-space where the ratios of the different invariant densities (taken 
pairwise) remain constant (equal to the respective initial values). 
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We note that because < kx^ + 5k^^ < 5k^, and because because the Fj0rtoft 
argument operates only with the strong inequahties and <^ rather than > and 



>), we can replace the zonostrophy density ( 14 1 with a simpler expression, 



Ck-^. (15) 

Thus we have for the boundaries between the different cascade sectors, see 
Fig. 

• E/Q boundary: As for the 2D Navier-Stokes turbulence considered before, we 
have for this boundary separating the energy and the enstrophy cascades, 

k^-ki (16) 

(i.e. a circle in the 2D A;-space, k^ + k^, — kg.) 

• E/Z boundary: Equating the ratio of the energy density \cOk\ to the zonostropy 
density kj' /k^ to the initial value of this ratio, we get for the boundary separating 
the energy and the zonostrophy cascades, 

k^/kx^kl/kox. (17) 

• Q/Z boundary: Equating the ratio of the enstrophy density k^ to the zonostropy 
density kx^ /k^ to the initial value of this ratio, we get for the boundary separating 
the enstrophy and the zonostrophy cascades, 

k^/kx ~ kl/kQx. (18) 

The first of these expressions, ( [T6] l, says that (like in 2D turbulence before) the 
energy must go to larger scales and the enstrophy must go to smaller scales. 



The second expression, (17i, says that the energy must go to the zonal scales, 
ky ^ kx- Moreover, this expression also poses a particular restriction on the path of 
the energy to the zonal scales, e.g. for ky ^ kx it should zonate at least as fast as 

ky = const k^J^ , see Fig. 

The last relation, ( [TS] ), is also interesting. Because k > kx, the curve ( 17 1 intersects 
the kx-axis at a finite distance, 

k:-kf/kli\ (19) 

We see that the zonostrophy cannot cascade too far to large ^'s, unless one starts with 
nearly zonal scales, koy ^ kox- In particular, if koy = kox we have k* = l^^^ko, i.e. the 
maximal allowed wavenumber for the zonostrohy cascade is practically the same as 
the initial scale. In other words, in this case the zonostrophy can only cascade to the 
larger scales. 
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5.3 Alternative argument for zonation. 



The generalized Fj0rtoft argument presented above was based on the zonostrophy 
conservation which was proven for the wave kinetic equation and therefore, it is 
expected to work well for the weakly nonlinear case. However, later we are going to 
present numerics in which we "push" the formal boundaries of validity and test the 
performance of the triple cascade picture in the case when the initial nonlinearity is 
strong. 

Because we are going to consider cases with strong nonlinearity, we would like 
to mention an alternative argument for zonation which makes sense when the ini- 
tial turbulence is strong. Suppose that turbulence is forced strongly at large A:'s, 
so strongly that the linear term is negligible compared to the nonlinear one in the 
CHM equation ([T]i. Since the linear term is the only source of anisotropy in the 
CHM model, the system is expected to build an isotropic inverse energy cascade 
identical to one of the 2D Navier-Stokes turbulence. As the inverse cascade pro- 
gresses toward the larger scales, the linear dynamics (j3-effect) become more and 
more important. A well-known boundary exists which defines the crossover from 
strong to weak turbulence, dominated by weakly nonlinear waves. This boundary is 
defined by the transitional wavenumber where the characteristic times of the linear 
and the nonlinear dynamics become equal. As it follows from the above speculation 
about the inverse cascade, the scaling for the nonlinear time has to be taken from a 
Kolmogorov-type estimate, X^l ~ (efe^)^'''^ and the linear time is just the inverse 
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wave frequency, Tl — k^/{j5kx). Equating Xnl and Tl we get 

k^=k^pkl, (20) 

where k^ — ()3^/e)'/^. This curve is plotted in Fig. (jsj) and it is the well-known 
'dumb-bell' or iazy 8' curve |'32lj33'|. For the modes which lie outside the curve, 
isotropic turbulence is the dominant process, whereas for those modes inside the 
dumb-bell, anisotropic Rossby wave turbulence is the dominant process. 




The alternative zonation argument is based on the observation that weak wave 
turbulence is much less efficient in supporting the energy cascade than the strongly 
nonlinear interactions. Thus one can suppose that the energy cascade does not pen- 
etrate inside of the dumb-bell curve, but rather, it turns and continues along this 
curve to the larger scales. This condition that the linear and the nonlinear times are 
balanced scale-by-scale could be called a "critical balance', in analogy with similar 
arguments about the critical balance state discussed in MHD turbulence 1,34] . 

The critical balance picture would mean that the energy cascade path of zonation 
would be given by expression (20 1, which incidentally is very close to the path 
predicted by the Fj0rtoft argument ( 17 1, if k^ ^ kp (the power 3/8 indeed very close 
to 1/3). 

A closer look reveals that the alternative zonation argument is actually not so 
"alternative" as it does not contradict our generalized Fj0rtoft argument based on 
the zonostrophy concentration. Indeed, if the energy followed a path through the 
scales at which the linear and the nonlinear timescales balance then the zonostrophy 
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would flow below this path, i.e. to the weakly nonlinear scales. Thus, the zonos- 
trophy would be supported at the weakly nonlinear scales and therefore would be 
conserved (even though both the energy and the enstrophy would be supported in 
the intermediately and strongly nonlinear scales in this case). But then all the three 
invariants are conserved (no need for weak nonlinearity for conservation of E and 
^2) and the generalized Fj0rtoft argument is once again valid. 

One should also keep in mind that for the dumb-bell curve the estimation of the 
shape is only approximate because it is based on the expression for T;vl which is 
valid, strictly speaking, only at fc's which are much greater than the dumb-bell (i.e. 
in the region of the isotropic inverse cascade). With this observation, we can put 
together the "kinematic" view of the Fj0rtoft argument and the "dynamic" critical 
balance condition, and suggest that: 

• Invariants E, Q. and Z cascade in the sectors prescribed by the generalized 



Fj0rtoft argument, i.e. ( 16 1, ( 17 1 and ( 18 i 



• The energy cascade follows a path along the scales at which Xnl ^ '^l which lies 
somewhere in the sector ( [T7] l and which does not strictly follow the dumb-bell 
prescription ( pO] ) (because of the unprecise definition of the dumb-bell discussed 
above). 

• This picture is of course preceded by the usual isotropic inverse cascade of E in 
the case when the initial turbulence is very strong. 

We emphasize that in this section we only considered the case when the initial 
turbulence is strong. If the initial turbulence is weak, i.e. well within the dumb-bell, 
then the energy cascade path remains in the weakly nonlinear scales (at least for a 
while) rather than lie on the dumb-bell curve. 

Now we will proceed to testing our theoretical predictions about the conservation 
of the zonostrophy and the triple cascade behavior via numerical simulations of the 
CHM equation in the Fourier space, Eqn (|2]i. 



6 NUMERICAL STUDY 

A pseudospectral code has been written to solve equation Q. No dissipation is used 
and the initial condition is given by 

rk|r=o-Ae^^+*'+ image, (21) 

where and k^, are constants and 0k are random independent phases, and by "im- 
age" we mean the mirror-reflected spectrum with respect to the A;, axis. Note that 
only the semi-plane > was used in our computations because of the symmetry 
i/Z-k = V^k arising from the fact that the streamfunction y^f in the CHM equation ([Tji 
is a real function. 

We opted to simulate such an evolving non-dissipative system rather than a 
forced/dissipated steady-state turbulence considered by the Fj0rtoft argument be- 
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cause it appears to be more physically relevant (since there appears to be no phys- 
ically meaningful dissipation acting selectively on nearly zonal and nearly merid- 
ional scales only). Yet, we hope that the cascade picture obtained via the Fj0rtoft 
argument is meaningful for such decaying turbulence too, similar to what appears 
to be the case e.g. in the 2D Navier-Stokes turbulence. 

Of course, when calculating a non-dissipative system one has to be aware of 
the possible bottleneck accumulation of turbulence near the maximum wavenum- 
ber after the turbulent front reaches these scales. Thus we make sure to stop our 
simulations before this happens. 



6.1 Centroids 

To quantify the cascades of the energy, enstrophy and zonostrophy in the time- 
evolving non-dissipative turbulence we introduce the centroids ("centres of mass") 
of the total of each invariant defined respectively as follows. 



Of course, it is not apriori clear if the Fj0rtoft argument formulated for the 
steady-state forced/dissipated turbulence has a predictive power for the trajectories 
of the centroids in the ^-space. In the present paper we "experimentally" verify that 
this is indeed the case. 

For the 2D Navier-Stokes (Euler), it is actually possible to recast the Fj0rtoft ar- 
gument for the non-dissipative evolving turbulence directly in terms of the centroids 
in a rigorous way, see Appendix B. However, the structure of the CHM is more in- 
volved and it is not clear if one can produce a generalized Fj0rtoft's argument for 
the triple cascades in terms of the centroids in a rigorous way. This is certainly an 
interesting question to be addressed in future. N 




(22) 




(23) 




(24) 



We have been able to find some of the relevant inequalities in terms of the centroids, but most of 
the conditions restricting the triple cascade sectors are still missing. 
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6.2 Settings for the weakly nonlinear and the strongly nonlinear 
runs. 



We have chosen two sets of parameters to be used in two runs corresponding to 
weak and strong initial nonlinearities respectively. 

The weakly nonlinear and the strongly nonlinear runs were performed at res- 
olutions 512^ and 1024^ respectively. This is because the weakly nonlinear sys- 
tems evolve much slower that the strongly nonlinear ones and one has to compute 
them for much longer Correspondingly, the centre of the initial spectrum and its 
width were chosen to be ko = (20,20) and A:* = 8 in the weakly nonlinear run and 
ko = (40,40) and k^,~ 16 in the strongly nonlinear run. 

The initial amplitudes in the weakly nonlinear and the strongly nonlinear runs 
are A — 10^^ and A = 5 x 10^^ respectively. Determining the relevant degree of 
nonlinearity a that corresponds to these initial conditions is tricky. One can directly 
estimate the linear and the nonlinear terms in the Eqn (j2]i and take into account the 
fact that there will be statistical cancellations in the sum of the nonlinear term due 
to the random phases, i.e. schematically 



^ individual-term^ 



I individual-term I 



This way we get an estimate 



(25) 



which gives a ^ 0.09 for the weakly nonlinear run and a ^ 0.7 for the strongly 
nonlinear run. However, it is likely that in the strongly nonlinear case the phases 
will quickly become correlated to a certain degree. Evaluating the nonlinear term 
in the extreme case when the phases are totally coherent would give an extra factor 
of Vn ~ 2^Tik\ ^ 100, so we would have for the strongly nonlinear run a ^ 70, 
which obviously is an overestimate. From common sense, in a non-rigorous way, we 
believe that the true relevant value for a in this case is closer to the random-phase 
estimate, but perhaps slightly higher, e.g. a ^ 1—2, which is approximately "on the 
dumb-bell". 



6.3 Weakly Nonlinear Case 

For this case, the degree of nonlinearity is weak, a ^ 0.09, so that the initial tur- 
bulence is well within the dumb-bell. Fig. Q shows the conservation of each in- 
variant. Because of the slow weakly nonlinear evolution, any quantity proportional 
to the turbulent intensity could look "conserved", so to demonstrate the true con- 
servation of the zonostrophy we plot its time evolution along with a non-conserved 
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quantity, / Xj/^dx. It is clear that all three invariants, E,Q and Z are well conserved. 
Namely, the energy is conserved within 0.01%, the enstrophy - within 0.15% and 
the zonostrophy is conserved within 1 %. Note that this is a first numerical demon- 
stration of the conservation of the zonostrophy invariant. We remind the reader that 
Z is precisely conserved by the wave kinetic equation (|9]l and therefore its conser- 
vation by the dynamical CHM equation (|2]) is subject to the applicability conditions 
of this kinetic equation, namely the weak nonlinearity and the random phases. It 
is not apriori clear how well these conditions are satisfied throughout the A;-space 
(particularly near the zonal scales). 

The cascade directions for £, i2 and Z are plotted in Fig. (jSj) in terms of the paths 
followed by the respective centroids, ( [22| , ( [23| and p4| . Note that for convenience 
we normalize the centroids to their initial values, A:^Q,fc^Q and fc^Q (which are differ- 
ent from ko) so that the centroid paths start from the same point in Fig. (jSj). 

In Fig. (j5]l we see that each invariant cascades well into its predicted sector In- 
terestingly, the enstrophy and the zonostrophy paths are well inside their respec- 
tive cascade sectors, whereas the energy follows the boundary of its sector with the 
zonostrophy sector. One should remember, however, that the boundaries between the 
sectors are not sharp because the Fj0rtoft argument operates with strong inequalities 
(^ and ^ rather than < and >). 




0.8 I ' ' ' ' 1 

10000 20000 30000 40000 
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Fig. 4 Conservation of energy, enstrophy and zonostrophy for the weakly nonlinear case. Non- 
conserved quantity f ^ also shown. 

Three successive frames of the energy spectrum in the 2D A:-space along with the 
x-space frames of vorticity distributions at the same moments of time for the weak 
nonlinearity case are shown in Fig. (j8]l. The initial spectrum, which represents a 
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Fig. 5 The cascades of energy, enstrophy and zonostrophy for the weakly nonlinear case, tracked 
by their centroids. 

gaussian spot centered at ko and its mirror image, is seen to grow "arms" toward the 
coordinate origin, so that a closed "ring" forms and then starts shrinking in size. The 
ring is suggestive of the dumb-bell shape (when complemented with the other half 
of the distribution at fc^ < 0), although the similarity is only visual rather than quanti- 
tative, because the nonlinearity is quite small. Presumably, the growing of arms and 
the ring shrinking are indicative of the structure of the anisotropic inverse energy 
cascade process. On the respective vorticity x-plots, we see initial dominant short- 
wave components propagating at ±45" (corresponding to the position of the initial 
maxima in the spectrum) which in time evolve into a more disordered turbulent state 
with a predominant zonal orientation. 



6.4 Strongly Nonlinear case 

For this case, a ~ 1 — 2 so that the initial turbulence is near the boundary of the 
dumb-bell. Fig. (j6]l shows the conservation of each invariant. While the energy and 
enstrophy are still well conserved (the energy within 0.2% and the enstrophy - within 
1.2%), the zonostrophy is not conserved initially. This is not surprising considering 
that the zonostrophy is only expected to be conserved if the nonlinearity is weak. 
What is more interesting, however, is that the zonostrophy growth saturates as time 
proceeds, so that the zonostrophy is rather well conserved in this case for large times. 
This suggests, as we argued before, that for large times the scales that support the 
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zonostrophy invariant are weakly nonlinear, even though the energy scales proba- 
bly remain moderately nonlinear, and the enstrophy scales are definitely strongly 
nonlinear. 

The cascade paths for £,12 and Z in terms of the respective centroids are plotted 
in Fig. (|7]i. Once again we see a picture which is similar to the one already observed 
for the weakly nonlinear case: the enstrophy and the zonostrophy cascades lie well 
inside their respective theoretically predicted sectors, and the energy cascade fol- 
lows the boundary of its sector. It is quite possible that the energy path lies in the 
"critically balanced" scales where the nonlinear and the linear time scales are of 
the same order However the measurement of the nonlinear time scale is quite am- 
biguous and it is still unclear if the critical balance approach can be formulated in 
a more precise way in this case. In any case, one can clearly see that even in this 
strongly nonlinear case the zonostrophy invariant is conserved for large times, and 
that the triple cascade picture predicted using this invariant provides a reasonable 
description of the turbulence evolution and explanation of the zonal jet formation. 

We have also looked at the energy spectra in the 2D wavenumber space and at the 
2D vorticity distributions on the x-plane evolving in time, see Fig. (|9]l. The essential 
features of the evolution of these distributions appeared to be remarkably similar 
to the ones of the weakly nonlinear case, with somewhat more evident zonal jets 
at later times. This can be explained by the fact that the strongly nonlinear systems 
evolve faster than the weakly nonlinear ones, so that what we see here is a more 
advanced stage of zonation. 
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Fig. 6 Conservation of energy, enstrophy and zonostrophy for the strongly nonlinear case. Non- 
conserved quantity f~ also shown. 
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Fig. 7 The cascades of energy, enstrophy and zonostrophy for the strongly nonlinear case, tracked 
by their centroids. 



7 SUMMARY 

In the present paper the generalized Fj0rtoft argument was used to predict a triple 
cascade behaviour of the CHM turbulence, in which the energy, the enstrophy and 
the zonostrophy are cascading into their respective non-intersecting sectors in the 
scale space. These cascades are anisotropic and the energy cascade is predicted to be 
directed to the zonal scales, which provides a physical explanation and the character 
of the formation of the zonal jets in such systems. 

The zonostrophy conservation, as well as the triple cascade picture, were tested 
numerically for the cases of both weak and strong initial nonlinearities. The zonos- 
trophy invariant was shown to be well conserved in the weakly nonlinear case. 
Moreover, the zonostrophy conservation was also observed for the case with strong 
initial nonlinearity after a transient non-conservative time interval. Presumably, this 
is because the zonostrophy moves in time to the scales that are weakly nonlinear 
(even though the energy and the enstrophy remain in the strongly nonlinear parts of 
the Fourier space). Using the energy, the enstrophy and the zonostrophy centroids 
for tracking the transfers of these invariants in the Fourier space, we demonstrated 
that all the three invariants cascade as prescribed by the triple cascade Fj0rtoft argu- 
ment in both the weakly nonlinear and in the strongly nonlinear cases. The energy 
appears to be somewhat special among the three invariants in that it tends to cas- 
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Fig. 8 Contour plots of the ID energy spectrum in A:-space (left) and corresponding contours of 
vorticity in x-space (right) for the weakly nonlinear case. 
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Fig. 9 Contour plots of the 2D energy spectrum in /:-space (left) and corresponding contours of 
vorticity in x-space (right) for the strongly nonlinear case. 



20 



Sergey Nazarenko, Brenda Quinn 



cade along the edge of the sector allowed by the Fj0rtoft argument, namely along 
the curve ^i^^. 

We beUeve that further studies would be helpful, both theoretical and numerical, 
for estabUshing the conditions under which the zonostrophy is conserved, in par- 
ticular finding out the extent to which the statistical properties of the system (e.g. 
random phases) are important in addition to weak nonlinearity of the zonostrophy 
supporting scales. It would be also interesting to study the behaviour of zonostro- 
phy in the other setups within the CHM model, e.g. the modulational instabiUty and 
truncated systems of coupled resonant triads. 
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Appendix A - Small-scale limit of Zonostrophy 

The general expression for the zonostrophy density is given by fTl\: 

L = arctan — arctan — (26) 

pk^ pk^ 

Expanding this expression in the powers of 1 /p up to ninth order, we get 

C- -2^/3^0+273-^'^ 



pfc2 p^k^ 

-2V3-^^[k^^+6klkf. + lkt] 

-2^/3 ^ [^kt + 27fc>2 + 15^2^4 ^ ^6] 



-2V3^^g^ [9/t° + mk°k;, + I26k^k'^. 
+2SkX.+k^.] + 0{p-^'^). 



22 Sergey Nazarenko, Brenda Quinn 

We note that the Taylor expansion of the frequency (O^ = —^/{k^ + P^^) in the 
powers of 1 /p is 

We see that in the leading order ^ = (0^, i.e. in the small-scale limit ^ is propor- 
tional to the energy and not an independent invariant. 

Thus, to find the truly independent invariant in the small-scale limit we must 
subtract this "energy" part. To have a simpler expression which is p -independent in 

the leading order, we will also multiply the result by — . Re-defined this way, the 
zonostrophy invariant is 

which in the Umit p — > oo becomes the expression we were looking for, 

C=^^^V^- (28) 



Appendix B - Fj0rtoft argument in terms of the centroids. 

Let us consider an evolving hydrodynamic 2D turbulence in the absence of forcing 
and dissipation. Here we will re-formulate the Fj0rtoft argument in terms of the 
energy and enstrophy centroids in the k and I (i.e. scale) spaces. This formulation 
will be rigorous and quite useful for visualizing the directions of transfer of of the 
energy and the enstrophy. In contrast with the version of the Fj0rtoft argument give 
in the main text, this formulation is for a non-dissipative turbulence rather than a 
forced/dissipated system. 

The energy and the enstrophy A:-centroids are defined respectively as 

/•OO 

kE = / kEkdk/E, (29) 

poo 

kn = / k^Ekdk/a, (30) 

Jo 

and the energy and the enstrophy /-centroids as defined respectively as 

Ie= ( k'^Ekdk/E, (31) 
Jo 

la^ [ kEkdk/£2=kEE/n, (32) 

JO 
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where Ej^ is the ID energy spectrum (i.e. the energy density in \k\). 

Theorem 1. Assuming that the integrals defining E^Q,kE,kQjE o.nd Iq converge, 
the following inequalities hold. 



kE < ^QJe, 
ka > yi27£, 
kEkn > Q/E, 

Ie > VE/n, 

In < VeJQ, 
IeIq >E/Q. 



(33) 
(34) 
(35) 
(36) 
(37) 
(38) 



We are going to prove this theorem using Cauchy-Schwartz inequality]^ which 
states that 



1/2 



1/2 



r f{k)g{k)dk < r f{k)dk r g^{k)dk 

Jo Jo Jo 

for any functions f{k),g{k) G L^. We will only deal with positive functions, so the 
absolute value brackets may be omitted. Being in in our case means that all the 
relevant integrals converge, as suggested in the statement of the problem. 

First, let us consider integral / kE dk and apply Cauchy-Schwartz inequality as 
foUows, 



kEdk= / {kE^I^){E^/^)dk< 



k^Edk 



1/2 



Edk 



1/2 



which immediately yields p3| ) and ([37]). 

Second, let us consider i2 — j k^E dk and split it as. 



k^Edk= {k^l^E^I^){k^l^E^I^)dk<( k^Edk 



1/2 



kEdk 



1/2 



which immediately yields ( 35 1. Combining ( [35) l with ( 33 1 gives ( 34 1. 
Now, let us split / kE dk in a different way. 



rkEdk= r{ky^E"^){k-'i^E"^)dk< { ( 

Jo Jo \Jo 



k^Edk 



1/2 



k-^Edk 



1/2 



which immediately yields ( p8| ). Combining ( [38] l with ([37]) gives ( [36| . 

We see that according to inequalities ( [33] l and ( [34| , during the system's evolution 
the energy centroid kE{t) is bounded from above and the enstrophy centroid kQ{t) 

^ The suggestion to use the Cauchy-Schwartz inequality for reformulating the Fj0rtoft argument in 
terms of the centroids was made to us by Gregory Eyink during the INI workshop the proceedings 
of which are published in this book. 
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is bounded from below (both by the same wavenumber k = y^£2/E), as one would 
expect from Fj0rtoft argument. Further, inequality ( 35 1 means that if kE (?) happened 
to move to small k's then kQ (f ) must move to large k's, that is roughly, there cannot 
be inverse cascade of energy without a forward cascade of enstrophy. Note that there 
is no complimentary restriction which would oblige kE{t) to become small when 
kn {t) goes large, so the A;-centriod part of the Fj0rtoft argument is asymmetric, and 
one has to consider the /-centroids to make it symmetric. Indeed, in additions to 
conditions ( 36 1 and (j37]l which are similar to ( 33 1 and ( 34 1, we have inequality ( 38 1 
meaning that if Iq (t) happened to move to small /'s then /^(f ) must move to large Z's, 
i.e. any forward cascade of enstrophy must be accompanied by an inverse cascade 
of energy. 

Importantly, we do not always have kE ^ I /Ie and kn ^ I /Iq . Indeed, consider 
a state with spectrum Ei^ ^ k^^l^ for k^ < k < kf, (with kf, 3> ka) and Ei^ = outside 

1 /3 2/3 

of this range. It is easy to see that for this state kE ^ ^b^a and Ie ^ l/kj, (i.e. 
kE oo 1/Ie) and kn k^. 



